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Abstract 

Decoherence effects on the traditional N vs. M photon coherent control of a two-level system are 
investigated, with 1 vs. 3 used as a specific example. The problem reduces to that of a two-level 
system interacting with a single mode field, but with an effective Rabi frequency that depends upon 
the fundamental and third harmonic fields. The resultant analytic control solution is explored for 
a variety of parameters, with emphasis on the dependence of control on the relative phase of the 
lasers. The generalization to off-resonant cases is noted. 
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I. INTRODUCTION 



Coherent control of atomic and molecular dynamics using optical fields has attracted 
much attention, both theoretically and experimentally Q,0,y]. Thus far, most theoretical 
work has focused on the idealized case of isolated systems, where loss of quantum phase 
information due to decoherence, i.e. coupling to the environment, is ignored. Such effects 
are, however, crucial to control in realistic systems, since loss of phase information results 
in loss of control. For this reason efforts to understand control in external environments 4j[- 
B and to _ t e for the result decohere |e. g . , fl-Q) ar e of great inte.es, 

There exist a number of basic interference schemes l|] that embody the essence of coherent 
control. One is the N vs. M photon scenario where control results from interference between 
state excitation using N and M photons simultaneously. In this letter we provide an analytic 
solution for control in the two-level N vs. M photon control scenario in the presence of 
decoherence. For simplicity, we examine the 1 vs. 3 photon case, although the solutions 
obtained below apply equally well to the N vs. M photon case, with obvious changes in the 
input Rabi frequencies and relative laser phases. 

In 1 vs. 3 photon control [16| a continuous wave electromagnetic field composed of a 
superposition of a fundamental and third harmonic wave is incident on a system. By varying 
the relative phase and amplitude of the fundamental and the third harmonic one can alter the 
population of the state excited by the incident field. Clearly, decoherence can be expected 
to diminish the 1 vs. 3 photon induced interference, and hence the control over excitation. 
A, t hou gh extent tneoretiea, Q - Q and expend, Q - Q studies have been 
carried out on the 1 vs. 3 photon coherent control scenario, there has been no serious 
examination of the stability of this control scheme in an external environment, barring a 
derivation of a simple analytical expression for the autoionization of a two-level atomic 
system for weak laser intensities, using the rate approximation [18]. Amongst the various 
possible influences of an environment on a system we focus on the loss of phase coherence, 
that is, dephasing. Dephasing is expected to occur on a time scale more relevant to control, 
since the duration of control field can be on the order of a picosecond or less, wheras the 
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typical time scale for energy transfer is considerably longer 

In this paper we show that the 1 vs. 3 photon phase control scenario (which controls 
the population) in a two-level system, when coupled to an environment, reduces to the 
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analytically soluble monochromatic field case, but with an effective Rabi frequency that is 
determined by the relative phase and amplitudes of the two fields. Sample results for control 
as a function of relative laser phase in the presence of dephasing are then provided. The 
possiblity of solving the off-resonance case is also noted. 

II. 1 + 3 PHOTON CONTROL 
A. Formalism 

Consider a two-level bound system interacting with an continuous wave (CW) electro- 
magnetic field and assume that the energy levels undergo random Stark shifts without a 
change of state during collisions with an external bath, e.g., elastic collisions between atoms 
in a gas. The CW field E(t) is treated classically, and the ground and the excited energy 
eigenstates states, of energy E\ and E 2 are denoted |1) and |2), respectively. 

In general, the system density operator p obeys the Liouville equation, 

| = -i[i/ W ,,]-^. ID 

Here Hit) = H atom + H- m t, where the free atomic Hamiltonian term is 

H atom = E 1 \l)(l\ + E 2 \2)(2\ (2) 

and the atom-field interaction term within the dipole approximation is 

H int = -E(t)[(l\d\2}\l}(2\ + (2\d\l}\2}(l\] (3) 

with electric dipole operator d. The second term in Eq. ([TJ, 1Z, is a dissipative term that 
can have a variety of nonequivalent forms associated with various master equations. Below 
we assume simple exponential dephasing of the off-diagonal p^-. 

In the simplest 1 vs. 3 control scenario, a two-level system is subject to the linearly 
polarized laser field: 

E(t) = hsje^e-^f +£ h e luJht e-^ h +c.c], (4) 

where £j is the real time-independent amplitude and 4>j is the phase of the corresponding 
field, with j = h, f. Here the subscripts /, h denotes the fundamental and its third harmonic, 
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and "c.c." denotes the complex conjugate of the terms that precede it. The fields have 
frequencies Uf and Uh — 3uf, chosen so that the third-harmonic and the three fundamental 
photons are on resonance with the transition from the ground state |1) to the excited state 
1 2). In the standard scenario , control is obtained by changing the relative phase 

and amplitudes of two fields, which results in the alteration of the degree of interference 
between the two pathways to the excited state. 

Within the rotating-wave approximation, the slowly varying density-matrix elements of 
the states |1) and |2), an = pa, (i = 1,2) and a 21 = p 21 e 3l ^ ft+ ^^ obey the following set of 



equations: 



^ = -\m[{pf 2 ] £}/h + p l2 £ h e^/h)a 21 ] 
^ = lm[(^£ 3 f /h + ^ 2 S h e^/h)a 21 ], 

^ = -7^21 + fa$>£f/h + P2i£ h e-^ /h)(a n - a 22 ), (5) 



with 



(3) _ J- Plnpnmpm2 

Ml2 ~ (2h) 2 ^ (u nl - Uf){u) f - uj 2m ) ' 

Here 7 P is the dephasing rate, u> nm is the frequency difference between levels \n) and \m) 
and jj, nm = (n\d\m). The quantities p\ 2 and /% denote the one-photon matrix element for 
the harmonic field and the effective three-photon matrix element for the fundamental field 
for the |1) — > 1 2) transition. Below, we use p = p\ 2 and p^ = pi 2 \ omitting the subscripts 
for simplicity. The controllable relative phase is (f) = <ph — 30/. 

It is convenient to define the one- and three-photon Rabi frequencies by Qh — ^£h/^ 
and VLf = p^£j/h, given in terms of their amplitudes and phases, by Vlh = \£lh\e l9h and 
f2j = \Vtf\e ie ! . Note that, although p and p( 3 ' are real for a bound system, we derive all 
the equations under the assumption that they can be complex so that the analysis can be 
extended to complex matrix elements arising in transitions to the continuum. Since £h and 
£f are real and positive, Oh and Of are determined by p and p^ 3 ': 

e ldh = Ai/M, (7) 
e ie f = p {3) /\p {3) \. (8) 

To amalgamate these Rabi frequencies and the relative laser phase of 0, we define the 
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effective Rabi frequency fl e s' 

Q ea e ie = n h e i(t> + n f = \Q h e^ + Q f \e ie , (9) 
where fl e s is real and positive. Here f2 e g- and 9 are related to flh and Qf as 



^eff = Vl^l 2 + l fi /l 2 + 2|fi fc fi/|cOS$, (10) 

\Q f \ \Q f \ 

tan# = [sin(0 + h ) + sin 9 f ] /[cos (<j> + 6 h ) + \-f\cos9 f ], (11) 

| "ft | \^h\ 

where 

$ = <t> + 9 h -e f . (12) 

It is worth noting some features of fi, e g that are evident from Eq. fjlOj) . First, the total 
excitation probability obtained in lowest order perturbation theory for 1 vs. 3 photon phase 
control in a two- level system [r| is proportional to f2 e fj. Hence, Q e s can be used to predict 
the controlled population, and its dependence on 0, when the fields are weak. Further (see 
below), Q e g plays a major role in determining the transient behavior of the excited state 
population for any field intensity in the absence or presence of dephasing. Second, the 
interference term in Eq. (fTU|) can be controlled by varying $, that is, by manipulating the 
relative phase <p of the two fields. Since /i and are real in a bound system, possible values 
of 9 h — 9j are and ±7T. When 9 h -9 f =0, i.e. , fifi^ > 0, cos $ = cos d>. On the other hand, 
when 9h~ Of = ±7r, i.e. , /i/^ 3 - 1 < 0, cos$ = — cos0. Thus, opposite interference effects are 
observed depending on the signs of \i and jt/ 3 '. Third, \\£lh\ ~~ 1^/11 — ^cff — \&h\ + 1^/1 
so that maximal interference effects occur when \Qh\ = |^/|- If \&h\ 7^ the smallest 
Q c s is not zero, and thus complete destructive interference, that is, zero excitation from the 
ground to the excited state, does not occur. 

Rewriting Eq. (jSJ) in terms of Q e g and 9, 

<Mu : -Im[fi cfr e l %i] (13) 



dt 

Im[fi cff e 4 V 21 ], (14) 



dt 

^ = -7 P cr 2 i + % -p eS e- ie \{a xl - o 22 ). (15) 
and introducing u = 2 Re(<7i 2 e~ ie '), w = 2 Iu^cx^e - * 9 ), and w = ct 22 _ gives 

du/dt = —-f p u, (16) 



dv/dt = — J p V + fl e ftW, 

dw/dt = —Q eS v. 



(17) 
(18) 



The resultant equations are now of standard form|28j, but with f2 e g replacing the Rabi 
frequency of the single field case discussed in Ref [28]. Note that the longtime steady- 
state solution to Eqs. (fTBj) to (JTSj) is found by setting du/dt = dv/dt = dw/dt =0, giving 
u(t — > oo) = v(t — > oo) = u?(t — > oo) =0. This implies that, regardless of initial conditions 
and for sufficiently large time, pure dephasing leads to an equilibrium state with equal 
populations in the ground and the excited states and with no remaining coherence. 

Substituting Eq. (jT7Jl into Eq. (jTHJl gives a simple equation for w : 

d 2 w /dt 2 + -f p dw /dt + ttl s w = 0. (19) 

In the important case where initially the ground state is populated and the coherence is zero 
[i.e., w(0) = — 1, u(0) = v(0) = 0], the excited state population p 22 = cr 22 is given by 

P22 = -^-[cos(si) + J sin(st)] + - for 7p < 20 eff , (20) 

<» - '" A ^!y'' - 5 f » % > < 21 ) 

e _2|£ . j 



where s = |y4£7p ff — 7^, and Ai j2 = §[— 7 P ± y7 p — 4f2g ff ]. The general behavior of the 
solution is seen to be determined by relative size of the dephasing time and the period of the 
Rabi oscillation. Analogous analytic results can be obtained for a 12 which decays with rate 
7 P . If the external field is intense enough so that 7 P < 2f2 e fj , then p 22 shows oscillations that 
are exponentially damped with time. On the other hand, if dephasing dominates over the 
Rabi oscillation, so that 7 P > 2Q c g or 7 P = 2f2 e ff, p 22 increases monotonically. However, in 
all cases p 22 reaches a stationary value of 0.5 at long times and p 22 ~ e | for short times. 



B. Sample Computations 

The behavior of the excited state population for several values of 7 P for a given value of 
QcS (here chosen as 2ir) is sketched in Fig. ^ For 7 P < 2Q e s, the introduction of dephasing 
increases the period of the oscillation and causes the amplitudes to decay as e 2 . Although 
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this is a CW laser field case, we can extract the result for the field being switched off at a 
specific time, i.e., a square pulsed laser which is on from t — to t — tf, by examining the 
population at time tf. (This assumes that there are no additional energy levels excited by 
the frequency breadth of the truncated CW source). Significantly, one can end up with an 
increased p 22 even for a larger dephasing, depending on the pulse duration. For example, 
assume that we turn off the field at t=l. If there is no dephasing, then the excited state 
population at t — 1 is and is thus less than that of any of the other cases with dephasing. 
On the other hand, for 7 p > 2f2 eff (here j p > 47r), there is no oscillation; p 22 just increases 
monotonically towards 0.5, where the system reaches the steady-state slower with increasing 
dephasing. If we were to consider a pulse rather than a CW laser field for this relatively 
strong dephasing case, the excited state population would be expected to increase up to 0.5 
with the increase in the pulse duration. 

Typical behavior of p 22 and of the 1 vs. 3 photon phase control profile (i.e., p 22 as a 
function of generic phase control variable $) for several values of Rabi frequencies and 7 P 
are shown in Fig. El Here we assume that the fields are abruptly turned off at the times 
indicated in the figure captions to produce a square pulse and the intensities are chosen so 
that \Qh\ — 1^/1) to enhance the interference effects. The effective Rabi frequency is then 
Q c g = \Qh\ ^2(1 + cos $). While $ = leads to a complete constructive interference of 
the two transition amplitudes, $ = 7r leads to a complete destructive interference, i.e. , no 
excitation from the ground to the excited state. 

The typical control behavior seen in Fig. El depends upon the pulse duration, as well as 
upon fl e g and j p . For example, when the field is weak and 7 P = 0, then p 22 is given by [Eq 

©-(EH)] 

P22 = T^t 1 _ cos(f2 cfr t)] for 7 P = (23) 
ss -|ftk| 2 (l + cos$)i 2 ; for Ismail. (24) 

Hence, the system shows a "cos $ rule" . [ A similar rule obtains from Eq. ([20"1) - Eq. ([22"ll 
when 7 P < 2fi cff and st « 1, when 7 P > 2fi cff and Ait and A 2 t are much less than one, and 
for 7 P = 2£7 eff when j p t/2 « 1.] Note also that Eq. (j22J) predicts oscillatory behavior of 
P22 as a function of Q e g at fixed t, as observed later below. 

The control profiles for small Qh and 7 P = (thin dashed lines in Fig. EJ) are then seen 
to be monotonically decreasing from the maximum excitation at $ = 0, to zero excitation 
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at <3? = it, i.e. they follow the "cos$ rule". By contrast, for strong intensity (thin solid 
lines in Fig. |2J) in which there are many Rabi cycles during the pulse, the control curve is 
not necessarily monotonic since the final excited populations are determined by the time at 
which the fields are turned off. Introducing dephasing is seen to lead to a decreased range 
of control whose magnitude depends on the relative strength of the dephasing and on the 
effective Rabi frequencies, according to Eqs. PO]) - 

Figure El demonstrates that phase control profiles are strongly dependent on the pulse 
duration. For weak intensities, as the pulse duration increases, the degree of control improves 
and the control curve continues to approximately follow a cos$ law (e.g., Eq. (|24jl ). This 
behavior is seen both in the absence and in the presence of dephasing, although dephasing 
reduces the yield for a given pulse duration. In the strong field case the control profile varies 
strongly with pulse duration. In particular, with j p — 0, if the pulse duration is smaller than 
the oscillation period (=1/2) of the fi cff ($ = 0) case, then p 22 decreases with increasing $, 
as shown in Fig. Efa). For the pulse duration greater than that period, the control profile 
no longer follows cos*!' and the maximal yields start to appear at $ 7^ 0, as shown in Figs. 
Efb) to (d). In all strong intensity cases, the addition of dephasing results in a decay of p 22 
with a rate of e ~ for a given fi e g. Thus the degree of the control worsens in the presence 
of dephasing as the pulse duration increases. Note that the introduction of dephasing leads 
to a degree of control C that converges to 0.5, where C is defined as the difference between 
the maximum and minimum excited state populations. 

Finally, we note that this treatment can be extended in two directions. The most obvious 
is to extend it to the general two-level iV photon + M photon interference scenario where 
the structure of the problem is exactly the same as that of the 1 vs. 3 photon case. The 
equations above therefore hold, but with the one and three photon Rabi frequencies and 
phases replaced by the N and M photon Rabi frequencies and phases. The second is to 
consider the more general case that includes the equal detuning of both fields from the 
|1) — ► |2) transition, i.e. 5 = uj 2 i — — ^21 — 3w/, and where the populations of levels |2) 
and |1) decay with the same rate jj. Then Eqs. ((Hj) to (fT5|) become: 



dan 
dt 
da 22 



7d(°"n - °"ie) - lm[Q eS e l9 cr 21 ] 



7d(tf22 - 02e) + lm[Q e{i e l a 2 i] 



(25) 



(26) 



dt 

9(721 

dt 



(id + lp + iS)a 2 i + -[ficfre 1 ](au - cr 22 ) 



(27) 



S 



Here o\ e and o"2 e are the steady-state values of o\\ and 022, respectively, when Q e g = 0, and 
are introduced to allow for relaxation to equilibrium. In terms of u, v, and w, the above 
equations lead to 

u 

du/dt = -Sv - — , (28) 
dv/dt = 5u- — + fiefrw, (29) 

dw/dt = - (w ~ We) + Q eS v. (30) 

where w e = cr 2e — ai e , T\ = 1/jd an d ^2 = l/(7d + 7p)- Note that these are then of the 
same form as the usual Bloch equations for a monochromatic field. Torrev gave detailed 
analytical solutions for these equations in the monochromatic field case |28l |29j and the 
same analytical solutions for the 1 vs. 3 photon phase control case can be used, where the 
single field Q considered by Torrey is replaced by Q c g. We do not pursue this direction in 
this letter. 



III. SUMMARY 



In summary, we have obtained an analytic solution for iV vs. M photon phase control of 
a two- level system in an environment described by a l/j p dephasing time, with N — 1 and 
M = 3 as a specific example. The results should serve as a prototype for understanding the 
results of N vs. M photon phase control in more complicated systems, such as controlled 
Xenon ionization and IBr photodissociationj^]. 
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Figure 1. Excited state population as a function of time for various dephasing rates, 7 P 
shown inside the box for Q e g = 2n. All the variables are in dimensionless units. Data points 
are connected by straight lines as a guide. 

Figure 2. Excited state population versus relative phase for two different Solid lines 
and dashed lines denote the case at = 2tt and \Qh\ — 7r/5, respectively. Thin lines and 
thick lines denote the case at 7 P = and 7 P = ir, respectively. Fields are turned off at (a) 
t = 0.25, (b) t = 0.25 x 2, (c) t = 0.25 x 3, and (d) t = 0.25 x 8. Note that the data in 
panel (d) is too widely spaced to produce the last P22 — 1 maxima, whose exact locations 
can be predicted from Eq. (|23J). All the variables are in dimensionless units. Data points 
are connected by straight lines as a guide. 
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